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Why polymer chains in a melt are not random walks 
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Abstract. - A cornerstone of modern polymer physics is the 'Flory ideality hypothesis' which 
states that a chain in a polymer melt adopts 'ideal' random-walk-like conformations. Here we 
revisit theoretically and numerically this pivotal assumption and demonstrate that there are no- 
ticeable deviations from ideality. The deviations come from the interplay of chain connectivity and 
the incompressibility of the melt, leading to an effective repulsion between chain segments of all 
sizes s. The amplitude of this repulsion increases with decreasing s where chain segments become 
more and more swollen. We illustrate this swelling by an analysis of the form factor F{q), i.e. 
the scattered intensity at wavevector q resulting from intramolecular interferences of a chain. A 
'Kratky plot' of q^F{q) vs. q does not exhibit the plateau for intermediate wavevectors characteris- 



tic of ideal chains. One rather finds a conspicuous depression of the plateau, 5{F ^{q)) 
which increases with q and only depends on the monomer density p. 
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Polymer melts are dense disordered systems consisting 
6i macromolecular chains. Theories that predict proper- 
ties of chains in a melt or concentrated solutions gener- 
cilly start from the 'Flory ideality hypothesis' [1]. The 
hypothesis states that polymer conformations correspond 
to those of 'ideal' random walks on length scales much 
larger than the monomer diameter [1-3]. The commonly 
Accepted justification is that mtrochain and mterchain 
excluded volume forces compensate each other in dense 
Systems [2,3]. This compensation has several important 
consequences. For instance, the radius of gyration R{s) 
of chain segments of curvilinear length s < N scales as 
i?^(s) = a^s, where N denotes the number of monomers 
per chain and a the statistical segment length of the chain. 
This result holds provided s is sufficiently large for all 
Ibcal correlations to be neglected. For s = A, it im- 
plies that the radius of gyration of the total chain obeys 
Rg = R{s = N) = a^/N [3]. A further consequence of 
chain ideality is that the intrachain scattering function, 

the 'form factor' F{q), is given by the Debye formula, 
= 2A(exp(-(gi?g)2)-l + (gi?g)2)/(gi?g)4 [3]. 

For intermediate wavevectors, l/i?g ^ q <C 1/a, F'^^\q) 
reduces to the power law F(o)(g) « 2/{qaf. Due to 
this power-law behavior, we refer to the latter g-regime 
as 'scale-free regime' in the following. 

Neutron scattering experiments have been extensively 
used to test Flory's hypothesis [4]. This technique allows 
one to extract F{q) from the total scattered intensity of a 
mixture of deuterated and hydrogenated polymers. One 
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Fig. 1: Role of incompressibility and chain connectivity in 
dense polymer solutions, (a) Sketch of the segmental corre- 
lation hole of a marked chain segment of curvilinear length s. 
Density fluctuations of chain segments are correlated, since the 
total density (dashed line) cannot fluctuate in dense polymer 
solutions, (b) Self-similar pattern of nested segmental correla- 
tion holes aligned along the backbone of a reference chain. The 
large dashed circle represents the classical correlation hole of 
the total chain (s = N) [2] . Here we argue that incompressibil- 
ity and chain connectivity lead to a repulsion of the segmental 
correlation holes, which increases with decreasing s. The en- 
suing swelling of chain segments manifests itself by deviations 
of the form factor from Debye's formula. 



usually plots q^F{q) versus q ('Kratky plot') to reveal the 
existence of a 'Kratky plateau' in the scale-free regime. 
In applications, however, the Kratky plateau appears to 
be elusive [4] . Possible causes for deviations are effects of 
chain stiffness and finite chain thickness. In special cases, 
the scattering signals from chain stiffness and thickness 
may compensate one another, leading fortuitously to an 
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extended Kratky plateau [5,6]. Apparently, 'Kratky plots 
have to be interpreted with care' [4]. 

It is generally believed that if all possible obscuring fac- 
tors are avoided and 'thin' flexible polymers are exam- 
ined, no deviations from the Kratky plateau should oc- 
cur. Here we show that this is not true because long-range 
correlations along the chain backbone are induced by re- 
pulsive interactions of chain segments in dense polymer 
systems [7,8]. This effect is related to the well-known cor- 
relation hole [2]. For chain segments of length 1 <C s < iV, 
the correlation hole leads to non-Gaussian deviations of 
the form factor. Our key claim is that for asymptotically 
long chains {N — > oo) the deviations are given in the scale- 
free regime by 
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with p being the monomer number density. In the fol- 
lowing, we first outline the theory (sketched in fig. 1) and 
provide a scaling argument for eq. (1). Computer simula- 
tions of two well-studied coarse-grained polymer models, 
both described in detail in Ref. [9], are then utilized to 
carry out a critical test of our predictions (figs. 2, 3 and 
4). 

The general physical idea of our theory and the long- 
range intrachain correlations it predicts is related to the 
'segmental correlation hole' of a chain segment of length s 
in d dimensions (fig. 1). Polymer melts are essentially 
incompressible (on length scales large compared to the 
monomer diameter), and the density p of all monomers 
does not fluctuate. On the other hand, composition fluc- 
tuations of marked chains or subchains may occur, sub- 
ject to the constraint of constant density. Figure 1 shows 
a marked chain segment of curvilinear length s within 
a much larger chain of length N ^ s. The monomer 
density distribution c(r, s) of the segment (dash-dotted 
line) becomes c* = c{r sa 0,s) w s/R{s)'^ close to the 
center of mass (r = 0) and decays rapidly at distances 
of order R{s). Since composition fluctuations are cou- 
pled by the density constraint, a second chain segment 
(of the same or another chain) feels an entropic penalty 
U*{s) w c*{s)/p when both correlation holes approach 
each other [7]. (The scaling of U*{s) is demonstrated be- 
low, fig. 2.) This repulsion when considered between two 
half-segments swells the segment. In d = 3, the effect is 
weak, U*{s) = s/pR{s)^ ^ ^/Vs, and a standard per- 
turbation calculation can be performed [3, 10]. This cal- 
culation considers quantities which are defined such that 
they vanish if the perturbation potential U* is switched 
off, and are then shown to depend, in leading order, lin- 
early on U* (which may be directly checked numerically) . 
For instance, one finds that [8] 
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describes the swelling of chain segments in real space. 
Note that a is the renormalized statistical segment char- 
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Fig. 2: Pair-distribution function p(r, s) of the centers of mass 
of chain segments of length s. Data for different chain and seg- 
ment lengths are successfully scaled by plotting u(r, s)/f/*(s) 
as a function of r/R{s) with U*{s) = s/pR{s)^. The figure 
compares different s for BFM chains of length A'' = 4096 and 
one BSM data set for A*' — 1024 (stars). The function is pos- 
itive for small r/R{s) where segments repel each other. It is 
weakly attractive for larger distances. For comparison we also 
indicate the self-density c{r, s)R{s)^ / s for one BFM example 
with s = 256 (dashed line) which is similar — albeit not iden- 
tical — to u(r^ s)/U*{s). 



acterizing infinite chains in the melt. In reciprocal space, 
this swelling appears in the form factor. As the segment 
length is related to the wavevector q and the ideal form 
factor F(°)(q) by s{q) - l/\q\^ - F^°\ we find for the 
deviations of the measured and ideal form factors 
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which agrees with eq. (1). The important feature of this 
|g|'^-correction is that it depends neither on the strength 
of the excluded volume interaction nor on the statistical 
segment length. Hence, it must be generally valid, even 
for semidilute solutions. We checked explicitly that the re- 
sult from the renormalization group theory for semidilute 
solutions [11, 12] takes the same form as eq. (1) with the 
amplitude 0.03124 . . . that is within 0.03% of our 1/32 [10]. 

The scaling of the effective interaction between 
two chain segments can be tested numerically by 
computing the correlation hole function u(r, s) = 
— ln(p(r, s)/p(oo, s)) where p{r^s) denotes the pair- 
distribution function of the centers of mass of all segments 
of length s. Since the correlation hole is shallow for large s, 
expansion leads to u{r, s) ~ 1 — p{r, s)/p{oo, s) « c(r, s)/p 
[7]. Hence the interaction strength at r/R{s) <C 1 should 
be given by U*{s) ~ u{0,s) ~ c*{s)/p. Figure 2 shows 
that a scaling in terms of the measured radius of gyration 
R{s) and U*{s) is indeed possible. This finding confirms 
that segments repel each other for small r/R[s) with an 
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Fig. 3: Kratky representation of the form factor F{q){qaY /2 
vs. rescaled wavevector qa. The axes are rescaled by the sta- 
tistical segment length a for the BFM and BSM (as indicated 
in the figure) to make both data sets readily comparable. The 
values of a for asymptotically long chains were obtained by a 
fit to eq. (2). Note that pa^ « 0.13 is similar for both models 
which explains the collapse for small q. (The agreement for 
large qa is fortuitous.) The Debye formula (dashed line), com- 
puted here for A'^ = 8192 using the measured radius of gyration 
(i?g ~ 120), overestimates the dip of the form factor at ga 1 
by about 20%. The bold line indicates the deviation from the 
Kratky plateau, eq. (1), with the choice F'°'(g) {qa)"^ /2 = 1, 
valid for ideal chains in the scale-free q-regime. The increase of 
F{q) for the largest wavevectors shown is due to non-universal 
correlations on the monomer scale. Such microstructure effects 
are not captured by the theory. 




Fig. 4: Scahng of the form factor deviations 5 {N/F{q)) /U* {N) 
plotted versus Q — qRg{N) as suggested by eq. (4). Here 
U*iN) = N/pR^iN) and S(l/F{q)) = 1/F{q) - l/F^^^q), 
where F(q) is the measured form factor and i^'"' denotes the 
Debye function. For each data set the scaling of the axes 
and the calculation of the Debye function are done with the 
measured radius of gyration Rg(N). Without free parameter 
perfect data collapse is obtained for all systems. Deviations 
from data collapse are expected for large Q due to finite per- 
sistence length and microstructure effects ('Bragg' regime). In 
the Guinier regime (Q ^ Qo), S{l/F{q)) increases rapidly as 
Q'^ (thin line). Most importantly, the non-analytic IQp-power 
law slope (bold line) can be seen over an order of magnitude 
in Q — a striking confirmation of the theory. Obviously, the 
larger the chains, the better the agreement with eq. (1) which 
is valid for asymptotically long chains only. 



entropic penalty U*{s) ~ I/a/s for large s ^ 10.^ 

Here, as in the subsequent figures, the body of our 
data comes from the bond fluctuation model (BFM) — 
a lattice model which we study by Monte Carlo (MC) 
simulations. Using a mixture of local, slithering snake, 
and double-bridging MC moves [9] we were able to equili- 
brate dense systems {p = 0.5/8) with chain lengths up to 
N = 8192 contained in periodic simulation boxes of linear 
size L — 256. Additionally, molecular dynamics simula- 
tions of a bead-spring model (BSM) [9] were performed to 
dispel concerns that our results are influenced by the un- 
derlying lattice structure of the BFM. For clarity, we show 
only one BSM data set {N = 1024, p « 0.84, L « 62). 
All length scales are given in units of the lattice constant 
or the bead diameter for the BFM or BSM, respectively. 
Compared to standard experimental polymers both mod- 
els have in common that the chains are very flexible and 
the units (beads) are point-like. They thus represent an 
ideal tool for testing the theory. 

Figure 3 presents a Kratky plot of F{q) obtained for 

^Strictly speaking, u{r,s) coincides with the effective interaction 
potential U* (r, s) of two independent segments of length s only if 
s = N. More generally, u{r,s) reflects also the interactions of their 
neighbors along the chains. However, this does not affect the scaling 
on short distances, r < R{s), which matters here. 



the longest chains currently available for both models. It 
reveals clear deviations from the ideal form factor (dashed 
line). Instead of reaching a plateau, the simulated F{q) 
continuously decreases with increasing q in the scale-free 
regime. This difference is expected from our theory. The 
result for inflnitely long chains — eq. (1) with l/F^^\q) — 
[qaY /2 (bold line) — enables us to understand the initial 
depression below the Kratky plateau. The prediction for 
infinite N cannot capture the decrease of the form factor 
for ga < 0.1 leading to the Guinier regime [4] where F{q) 
is determined by the finite size of the simulated chains. 
A clearer evidence for the theory should thus be obtained 
by a different comparison between theory and simulation, 
which accounts for finite- iV effects. 

This is achieved in fig. 4 which focuses on deviations 
5{l/F{q)) = 1/F - 1/F(o) from the Debye formula F^°\ 
If the axes are rescaled according to 

m{Q) = 6{N/Fiq))/U*{N) (4) 

with Q = qRg{N), a perfect collapse of the simulation 
data for all N and both polymer models is obtained by 
using only one length scale, the (independently) measured 
Rg{N). Eq. (4) is suggested by the fact that F{q) is pro- 
portional to the number of scatterers, TV, and is a func- 
tion of qRg{N) only [3]. The master curve m{Q) thus 
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generalizes eq. (1) to finite N. In the Guinier regime 
(Q <C Qo ~ 6.3) we find m{Q) ~ Q'*, as expected from a 
standard analytic expansion in Q^. (The first two terms 
in and must vanish by construction. ) The scale- 
free regime is observed for Q > Qo- Since S{l/F) ~ 
in this regime, it follows from eq. (4) that m{Q) ~ IQI"*, 
which provides an additional scaling argument for eq. (1). 
Fig. 4 demonstrates for the largest chains available over 
more than an order of magnitude in Q a good quantita- 
tive agreement (bold line). We therefore suggest that an 
analysis of experimental data should employ this scaling 
plot as a diagnostic for the effects discussed here. 

In summary, we have identified a general mechanism 
that gives rise to deviations from ideal chain behavior in 
polymer melts, even if the chains arc very flexible and thin. 
This mechanism rests upon the interplay of chain connec- 
tivity and the incompressibility of the system, which gen- 
erates an effective repulsion between chain segments of 
curvilinear length s (fig. 1). A polymer in dense solutions 
may not be viewed as one soft sphere (or ellipsoid) [13, 14], 
but as a hierarchy of nested segmental correlation holes 
of all sizes aligned and correlated along the chain back- 
bone. The repulsion scales like U*{s) ~ for s 1 
(fig. 2); it is strong for small s, but becomes weak for 
s N in the laxge-N limit. ^ The overall size of a long 
chain thus remains almost 'ideal', whereas subchains are 
swollen (eqs. (1,2)). This swelling manifests itself in the 
form factor. In the Kratky representation, F{q) displays 
a hump followed by a decrease in the (/-regime where the 
Debye function would exhibit a plateau (fig. 3). It should 
be stressed that eq. (1) does not predict an analytic ex- 
pansion in terms of as one might naively anticipate. 
The intriguing [(/^-asymptote for infinite chains (fig. 4) 
formally arises from dilation invariance of the calculated 
diagrams [10]. Established theoretical methods [15,16] as- 
suming implicitly analytical properties of scattering func- 
tions must therefore overlook these non-analytical terms. 

These deviations from ideality should be measurable by 
neutron scattering experiments of flexible polymers. Un- 
fortunately, finite persistence length effects may mask the 
predicted behavior if the chains are not sufficiently long. 
An experimental verification — e.g. following the lines 
of a promising recent study [6] — would be of great fun- 
damental interest; it could also delineate the conditions 
where the predicted corrections to ideality can be observed 
in real polymer systems and must be considered in under- 
standing their structure, phase behavior, and dynamics. 
Such an impact on structure and dynamics is expected 
theoretically, for the bulk [8, 17, 18] and also for thin poly- 
mer films [7]. Quantitative applications of classical the- 
ories, such as the Rouse or reptation models [3], may be 
plagued by systematic errors due to the neglect of these 

^Consequently, a contribution U*{N) ~ 1/Vn to the free en- 
ergy of a polymer chain in the melt exists scaling non-linearly with 
molecular weight. This may be demonstrated computationally from 
the size-distribution in systems of linear self-assembled equilibrium 
polymers. 



effective repulsive interactions. They could, for instance, 
be responsible for observed deviations from Rouse behav- 
ior [19]. Moreover, for thin polymer films the repulsive 
interactions are known to be stronger than in the bulk [7] 
providing, hence, a mechanism to explain the systematic 
trend towards chain swelling observed experimentally for 
film thicknesses below the bulk radius of gyration [20]. 
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